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We report the first evidence for the r)b{2S) using the hb(2P) — > r/b(2S)'y transition and the first 
observation of the ht(lP) — > rjt(lS)'y and ht(2P) —¥ rjb(lS)y transitions. The mass and width 
of the rib(lS) and rib(2S) are measured to be m Vb (i S ) = (9402.4 ± 1.5 ± 1.8) MeV/c 2 , m Vb ( 2 s) = 
(9999.0 ± 3.5± 2 j)MeV/c 2 and r %(is) = (10.8 Z.f° -to) MeV. We also update the h b (lP) and 
hb(2P) mass measurements. We use a 133.4 fb _1 data sample collected at energies near the T(5S) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 14.40. Pq, 13.25. Gv, 12.39.Pn 



Bottomonium, a bound system of a bb quark-antiquark 
pair, is well described by nonrclativistic quantum me- 
chanics due to the slow motion of the heavy quarks 
The spin-singlet (5" = 0) states of the quark pair with 
zero orbital momentum, L = 0, are customarily called 

r]h{nS) with n = 1, 2, The hyperfine splitting AA/hf 

from the corresponding spin-triplet state, i.e. the mass 
difference between the Y(nS) and rjh(nS), provides an 
important measure of the spin-spin interaction between 
the quark and the antiquark. The existing measure- 
ments Q of the r)h(lS) mass are only in marginal agree- 
ment with theoretical expectations [!, H| and the r/b(lS) 
width is yet to be measured. There is no information 
available on the radially excited state r/h(2S). 

In this Letter, we report the first evidence for the 
ijb{2S) in the hh(2P) —> rjb{2S)^ transition and the first 
observation of the h b (lP) — > i]b{lS)j and /i&(2P) — > 
r]b(lS)"f transitions. We use a 121.4 fb -1 data sample at 
the T(5S*) resonance and 12.0 fb _1 of energy-scan data 
collected nearby with the Belle detector [H[ at the KEKB 
asymmetric-energy e + e~ collider Q. 

We study the processes e + e~ — > T(55) — > 
hb(nP)-K + tt~ — > [r]b(mS')7]7r + 7r _ in which the rjb{mS) 



states are reconstructed inclusively. The approximate 
values of the expected energies of the photons in the 
hb(nP) rest frame are given in Table fl] The hb(nP) 
signal is tagged using the missing mass of the 7r + 7r~ 
pair A/ m i SS (7r + 7r~), while the r]b(mS) signal is tagged 
using the variable M^ ss (tt + 7r~ 7) = A/ miss (7r + 7r~7) - 
M m i ss (7r + 7r~) + m/, s („p). The missing mass is defined 
via M misB (X) = y/(E c . m . - E* x ) 2 -p*x, where E c . m . is 
the center-of-mass (cm.) energy and E\ andp^c are the 
energy and momentum of the system X measured in the 
cm. frame. We fit the M m i ss (7r + 7r~) spectra for different 

^miL( 7r+7r_ 7) bins to measure the hb(nP) yield. This 
procedure removes the background due to random 7r + 7r~ 
combinations. The hb(nP) yield peaks at M^l^ir+ir^j) 
values corresponding to the masses of the rjb(mS) states 
from the hb(nP) — ¥ rjb(mS)^/ transitions. 

The selection criteria for the 7r + 7r~ pairs are the same 
as those described in Ref. Q- We use events that pass 
the Belle-standard hadronic event selection and consider 
all positively identified tt + tt~ pairs that originate from 
the vicinity of the interaction point. Belle previously 
observed that the decay T(5S) — > hb(nP)ir + Tr~ pro- 
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TABLE I: Expected photon energies in h b {nP) — > rjbirnS)^ 
transitions (E-y), it veto parameters (AM and E t h) and the 
parameter a of the M^2. B (ir + tt~ 7) resolution function. 





S 7 , 
MeV 


AM, 
MeV/c 2 


MeV 


0, 

MeV/c 2 


h b (2P) - 




260 


10 


125 


13.6 ± 1.1 


h b (lP) - 


* 7ft(lS) 7 


500 


13 


75 


19.8 ± 1.1 


h b {2P) - 


* %(1S) 7 


860 


17 


75 


8.6 ±0.7 



ceeds via the intermediate resonances Zf,(10610) and 
Z & (10650) We exploit this with the additional re- 
quirement 10.59 GeV/c 2 < M iniBS (7r ± ) < 10.67 GcV/c 2 , 
which suppresses the combinatorial background by a fac- 
tor of 5 [1.6] for the h b (lP) [h b {2P)\ without any signifi- 
cant loss of the signal. Photon candidates are clusters in 
the electromagnetic calorimeter that are not associated 
with charged tracks. We apply a veto on 7r° — > 77 de- 
cays, rejecting a photon candidate if the invariant mass of 
it and any other photon in the event with energy above 
the threshold E t h is within AM of the ir° mass. The 
parameters -Eth and AM, listed in Tabic U arc chosen 
by maximizing the ratio where S is the number of 
signal events in the Monte Carlo (MC) simulation and 
B is the number of background events estimated from a 
small fraction (0.1%) of the data. To suppress contin- 
uum e + e~ — > qq (q — u, d, s, c) background, we use the 
ratio i?2 of the second- to zeroth-order Fox- Wolfram mo- 
ments 0. In the ?7b(15') analysis, we require R2 < 0.3, 
which was optimized using the T(5S) — > T(2S)tt + tt~ 
decays 0. For the high statistics /if, (IP) — > rjb(lS)~/ 
transition, the optimum is shifted to R 2 < 0.32, which 
we adopt here for the rjb(2S) analysis. 

To calibrate the photon energy resolution function, 
we use three control channels: D* a — > -fD°(^s- K~tt + ), 
7t° — > 77 and 77 — > 77. For the two-photon final states, we 
require that the energies of the photons in the laboratory 
frame be almost equal: \Ei -E 2 \/(E 1 + E 2 ) < 0.05. This 
reduces the resolution shape dependence to a single vari- 
able. The resolution shape is parameterized by a double- 
sided bifurcated Crystal Ball function [l(J in which a 
bifurcated Gaussian is smoothly joined with power law 
tails on both sides. The signal is extracted using the 
M(K-TT + ~{)-M{K-iT + ) + m D o distribution for the D*° 
and M(77) for the 7r° and rj. From comparisons of the 
peak positions and widths in data and MC simulation, we 
determine shifts in the photon energy AE/E and width 
correction factors / as a function of E. Various calibra- 
tion channels give consistent results; an uncertainty is 
assigned based on their spread. These translate to typi- 
cal mass shift and width-correction factors of 1.2 MeV/c 2 
and 1.13, respectively. Average values of left and right 
widths (tr) of the bifurcated Gaussian components of the 



We update previous Belle measurements of the hb(nP) 
masses [7[, incorporating a 10% increase in statistics and 
the requirement of an intermediate Zj, [ll|. The results 
of the fits are shown in Figs [T] and [5] and the fitted signal 
parameters are listed in Table HU The confidence level of 
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FIG. 1: (colored online) The M m i ss (7r + 7r~) spectrum in the 
hb(lP) region. In (a) the data are the points with error bars 
with the fit function (blue solid curve) and background (red 
dashed curve) overlaid, (b) shows the background subtracted 
data (points with error bars) while the signal component of 
the fit is overlaid (blue curve). The background is combina- 
torial. 




M miss (*V), GeV/r/ 



FIG. 2: (colored online) The M m i ss (7r + 7r~) spectrum in the 
ht(2P) region. The legend is the same as in Fig. [T] The 
background components are random combinations and K% — > 
tt + -k~ decays. The combinatorial component is shown in (a) 
by the black dotted curve. 



M^"g S (7r + 7r 7) resolution functions are given in Table HI the fit in the /ih(lP) [hb(2P)\ region is 35% [70%]. The es- 
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TABLE II: The yield and mass of peaking components from 
the fits to the M m i SB (-7r + 7r~). Here and everywhere in this 
Letter, the first quoted uncertainty is statistical (unless stated 
otherwise) and the second (if present) is systematic. 





N, 10 3 


Mass, MeV/c 2 


T(5S) - 
T(3S) - 
T(5S) - 


> h b (lP) 
■+ T(1S) 
■+ T(2S) 


70.3 ± 3.3±£;? 
13 ±7 
61.3 ±4.1 


9899.1 ±0.4 ±1.0 
9973.0 
10021.3 ±0.5 


T(55) - 
T(5S) - 
T(25) - 
T(5S) - 


■> T(1D) 
> h b (2P) 
■» T(15) 
+ T(35) 


14 ±7 
89.5 ± 6.11^ 
97 ± 12 
58 ±8 


10169 ±3 
10259.8 ± 0.5 ± 1.1 

10305.6 ± 1.2 

10357.7 ± 1.0 



timation of systematic uncertainties follows Ref. Q ■ For 

■%iK)> we find 



L Xb.j(nP) ■ 



hyperfine splittings ^ 2,7+ 1 i 
,/=o 

AMhf(IP) = (+0.8 ± l.l)MeV/c 2 and AM HF (2P) = 
(±0.5±1.2) MeV/c 2 , where statistical and systematic un- 
certainties in mass are added in quadrature. 

We fit the M m i SS (7r + 7r _ ) spectra for each 



^miss( 7r+7r 7) b m t° measure the hb(nP) yield as 
a function of M^ > ss (tt + 7t~j). We fix the masses of the 
peaking components at the values given in Table [TT1 
Although the M m i ss (7r + 7r _ ) combinatorial background 
shape in the hh(2P) region is rather complicated 
(described by an eighth-order Chcbyshcv polyno- 
mial), it changes slowly with M^ ga (ir + 7r~ 7) . For the 
M m i ss (TT + iT~) fits for each bin, we multiply the polyno- 
mial shape with parameters fixed at their values from 
the overall fit by a lower-order polynomial with floating 
coefficients; this procedure improves the accuracy of 
the hb(nP) yield measurements. For this lower-order 
function, we use a first- [fifth-] order polynomial for 

the rj b (lS) [776 (25)] Mi 2 i ) ss ( 7r+7r "7) region. From a 
generic MC simulation, we find that the K s — > tt + it~ 

contribution is independent of the ^mis S ( 7r+7r ~7) value 
in the ?7f,(15) region; in the 776 (25) region, we restrict the 
Mm SS (7T + fit range to 10.10 GcV/c 2 - 10.34 GcV/c 2 , 
thereby avoiding the sharp rise in the K s — > 7r + 7r~ 
contribution that occurs at 10.37 GeV/c 2 . The results 
for the hbilP) and hb(2P) yields as a function of 

-^miss( 7r+7r ~7) are presented in Fig. 03 Clear peaks at 
9.4GeV/c 2 and 10.0GcV/c 2 are identified as signals 
for the rj b (lS) and rjb(2S), respectively. Generic MC 
simulations indicate that no peaking backgrounds are 
expected in these spectra. 

We fit the hb(nP) yield dependence on M^ ss (n + 7r~ 7) 
to a sum of the rjb{nS) signal components described by 
the convolution of a non-relativistic Breit-Wigner func- 
tion with the calibrated resolution function described 
above and a background parameterized by an exponenti- 
ation of a first- [second-] order polynomial in the 776(16*) 




9.7 9.8 



9.9 10 10.1 



MftLfrVy), GeV/c 2 



FIG. 3: The h b (lP) yield vs. M^ ss (n + tt - 7) (a), and h b {2P) 
yield vs. M^ as (Tv + iT~'y) in the r] b (lS) region (b) and in the 
rjb(2S) region (c). The solid (dashed) histogram is the fit 
result (background component of the fit function). 



[r7&(25)] region. The two M^2, s (ir + ir~~f) spectra [from 
the hb(lP) and h b (2P)] with rjb{lS) signals are fitted si- 
multaneously. We take into account the effect of multiple 
photons in the same M^ ss (TT + n~"f) bin by increasing the 
errors in the corresponding Af m i SS (7r + 7r _ ) histogram. We 
find event yields for the hb(nP) — > r)b(mS) transitions of 
N 1P ^ 1S = (23.5±2.0)xl0 3 , N 2P ^is = (10.3±1.3) x 10 3 
and N 2P ^2s = (25.8 ± 4.9) x 10 3 ; the fitted masses 

(9402.4 ± 1.5 ± 1.8) MeV/c 2 , 
1%) MeV and m^ (2S ) = (9999.0 ± 
3.5l'i%)MeV/c 2 . The confidence level of the ^(15) 
[7/5 (25)] fit is 61% [36%]. If the r] b {2S) width is al- 
lowed to float in the fit, we find IV^si = (4l~2o)MeV 



»2P-*-2S = 

and width are m^ng) 

— I'm 0+4.0+4.5 
(is) - UU.s_ 3 7 

■ +2.8\ 



r 



t ^„^ 9t? ^ < 24McV at 90% C.L. using the Feldman 



776 (2S) 

Cousins approach [l2j|. For the ^(25) mass determi- 
nation and yield measurements quoted above, we fix 
the rjb(2S) width at its value from perturbative calcula- 
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tions T Vb{2S) = r i)6(ls) ^p} = (4.9l^)MeV, where 
the uncertainty is due to our experimental uncertainty in 
T n b (is)- _ 

To estimate the systematic uncertainties in the rjb(nS) 
parameters, we vary the polynomial orders and fit inter- 
vals in the M miss (7r + 7r~) & ^miL( 7r+7r ~7) nts > an d the 
^miss( 7r+7r ~7) binning by scanning the starting point of 
the 10MeV/c 2 bin with 1McV/c 2 steps. We also multi- 
ply the non-relativistic Brcit-Wigncr function by an P 3 
term expected for an electric dipole transition and in- 
clude the uncertainty in the /if, (IP) and hb(2P) masses 
and in the estimated value of the r)b(2S) width. The 
contribution of each source is given in Table IIIIl We 
add the various contributions in quadrature to estimate 
the total systematic uncertainty. We repeat the analy- 
sis using signal MC simulation instead of data and find 
no shift of the ijb(nS) parameters compared to the MC 
input. For the hyperfme splittings m X („s) — m ri b (nS) 
we determine AM KF (1S) = (57.9 ± 2.3) McV/c 2 and 
AM H f(25) = (24.3±4 °) MeV/c 2 , where statistical and 
systematic uncertainties in mass are added in quadra- 
ture. 



TABLE 111: Systematic uncertainties in the /ft (15) and 
ri b (2S) masses (in MeV/c 2 ) and in the r? 6 (lS) width (in MeV). 





m Vb (is) 


r %(is) 




M miss (7r+7r~) fits 


+0.1 

-0.2 


+0.0 
-0.8 


+0.0 
-0.1 


M i'iss( 7r+7r ~7) binning 


+0.2 

-0.5 


+0.7 
-0.8 


+2.3 

-1.4 


M^ii^-j) fit 


+0.2 
-0.0 


+4.2 
-0.6 


+0.5 
-0.2 


Calibration 


±1.4 


±1.5 


±0.7 


r)b(nS) line shape 


+0.6 

-0.0 


+0.0 
-0.4 


+0.6 
-0.0 


h b (nP) mass 


±0.9 


±0.0 


±1.1 


rj b (2S) width 






+0.1 

-0.0 


Total 


±1.8 


+4.5 

-2.0 


+2.8 

-1.9 



Using Wilks' theorem [15j, we find 15 a [9 a] for the 
h h (lP) -)■ 77b(lS)7 [h b (2P) -> 176(15)7] statistical signifi- 
cance. For the significance of the r)b(2S) signal, we use a 
method that takes into account the trial factor associated 
with the definition of the search window or the so-called 
"look-elsewhere effect" [HI]. To determine this window, 
we conservatively assume r = and r = 1 for the ra- 
tio r = Aj\/HF(is) • (F° r reference, the measured value 
of r = 0.4201q q7q is consistent with perturbative calcu- 

™ 2 . t^T(2S) 

0.513 ± 0.011, 



1 m 2 pT(2S) 

lations [14| that predict J (2S> 

m T(is) r ee 

where the error is due to the uncertainties in r ee .) We 
find the significance of the rjb(2S) signal to be 4.8 a (4.2 a 
including systematics) . 

Branching fractions B[ht(nP) — > -q^mS)^] are de- 
termined from K°^ mS /(N^ al e), where N?£ a \ 



N nP^mS K ti_CUt ] is th e number of the h b {nP) -> 
Tjb(mS)j transitions [hb(nP)] that are rejected by the R2 
and 7T° veto requirements [by the i?2 requirement], e is 
the reconstruction efficiency of the photon. In this way, 
we do not rely on MC simulation for the determination 
of the efficiency of the R2 and 7r° veto requirements. To 
determine N a f l ~ cut and N^Z^g, we repeat the analy- 
sis for the events rejected by the corresponding require- 
ments. In the M^ sa (ir + n~ 7) [M miss (7r + 7r~)] fits, we fix 
the r)b(mS) mass and width [the hb(nP) mass]. We es- 
timate the iV^p tI-cut systematic uncertainty by varying 
the polynomial order and fit range in the M m i ss (7r + 7r~) 
fit and by varying the hb(nP) mass within its uncer- 
tainty. We find 7Vto, tal = (84.2 ± 4.4 +H) x 10 3 and 



Ajtotal 
JV 2P 



(98.5 ±8.1 ±|;|) x 10 3 . For N tota \ we add the N 
and _/v antl_cut uncertainties in quadrature. The contribu- 
tions from the various sources of systematic uncertainty 
to the A^p^s are given in Table [TV] The h h {nP) signal 
shape uncertainty is not considered for the yields, since 
it cancels in the ratio N^pt\ mS / N^p &l . In combining 
the N nP ^ m s and ^pi^m's uncertainties, we take into 
account correlations of the calibration and rjb{nS) line 

(30.9 ±3.2 1 3 *) x 10 3 , 
16.1 ± 2.4±|;g) x 10 3 and N$g* as = (35.6 ± 



shape sources. We find N[p^ is 

Artotal 
JV 2P^1S 

7.3 It*) x 10 3 . 



TABLE IV: Systematic uncertainties in the total yields of 
the ft* (IP) -)• Vb(lS)j, h b (2P) -> i7 6 (15)7 and h b (2P) -> 
77b (25)7 transitions (in 10 s ). 





a rtotal 
"1P-*1S 


a rtotal 
JV 2P->1S 


Artotal 
JV 2P->2S 


M miss (-7T + n ) fits 


+0.2 
-1.1 


+ 1.6 
-1.8 


+0.0 
-1.3 


M i"L( 7r+7r ~7) binning 


+ 1.0 

-0.3 


+0.5 

-0.1 


+ 1.5 
-3.7 


Mil(7r + 7r- 7 ) fit 


+3.2 
-1.2 


+ 1.0 
-1.3 


+2.5 
-2.1 


Calibration 


+0.3 
-0.4 


±0.0 


±1.3 


rjb(nS) line shape 


+0.2 

-0.0 


+0.1 
-0.0 


+0.3 
-0.0 


r)b(nS) mass & width 


+0.7 

-0.5 


±0.3 


+3.1 
-2.8 


Total 


+3.4 
-1.8 


+2.0 
-2.2 


+4.5 
-5.4 



The efficiency is determined using phase-space MC 
simulated events that we weight according to the expec- 
tation 1 — cos 2 7 + 2 cos 6? 7 cos cos^ 7r + 7 [17[, where 



N, 



P-traS 



and N^ 1 



P->mS — 
anti— cut 

P ' 



Oy (#„+) is the angle between the beam axis and 7 (ir + ) 
momentum and 0^+^ is the angle between 7 and ir + mo- 
menta, with all momenta measured in the cm. frame. Ef- 
ficiencies to reconstruct the photon after the 7r + 7r~ pair is 
already reconstructed are 74.6% for /if, (IP) — > r/-f,(lS) 7 , 
73.4% for h b (2P) -> 176(15)7 and 76.1% for h b (2P) -> 
i~lb{2S)j. The r]b(lS) [rjb(2S)] signal function is normal- 
ized to unity in the mass window 9.3 GeV/c 2 — 9.5 GeV/c 2 
[9.9GeV/c 2 - 10.1GeV/c 2 ]; thus, the yields and effi- 
ciencies correspond to this mass window. The effi- 
ciencies have relative uncertainties of 2% due to pos- 
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siblc differences between the data and MC simulation, 
and ±08% [±06%] duc to tue uncertainty in the T) b (lS) 
[^(25)] width. Wc find B[h b (lP) —} t]b(lS)i\ = (49.2 ± 
5.7±|J)%, B[h b (2P) —} r) b (lS)i] = (22.3 ± 3.8 ±H)% 
and B[h b {2P) -> r] b (2S)-f} = (47.5 ± 10.5 These 
branching fractions are a factor of 1.2 to 2.5 higher than 
theoretical expectations [l8j . 

In summary, we report the first evidence for the rj b {2S) 
using the h b {2P) — > ?] b (2S)j transition, with a sig- 
nificance, including systematics, of 4.2 er, and the first 
observation of the h b (lP) — > T] b (lS)j and h b (2P) — > 
rj b (lS)"/ transitions. The mass and width parameters of 
the r) b (lS) and r] b (2S) are measured to be m^ng) = 
(9402.4 ± 1.5 ± 1.8) MeV/c 2 , m Vb[2S) = (9999.0 ± 
3.5± 2 J)MeV/c 2 and T m{is) = (10.8 ±f;° ± 4 2 .o) MeV. 
The m^ 6 (2S) an< i F^ng) are first measurements; the 
m^ns) measurement is more precise than the current 
world average and is (11.4 ± 3.6) MeV/c 2 above the cen- 
tral value [l9(. The hyperfine splittings, AMhf(IS') = 
(57.9 ± 2.3) MeV/c 2 , AM HF (2S') = (24.3±t°) MeV/c 2 
and their ratio 0.420±n nla , are in agreement with the- 
oretical calculations [4j, [20J. We update the h b (lP) and 
h b (2P) mass measurements mh b np) — (9899.1 ± 0.4 ± 

1.0) MeV/c 2 , m hb( 2P) = (10259.8±0.5±1.1) MeV/c 2 , and 
IP and 2P hyperfine splittings AM HF (1P) = (+0.8 ± 

1.1) MeV/c 2 , AM HF (2P) = (+0.5 ± 1.2) MeV/c 2 . These 
results supersede those in Rcf. pj . 
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